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Abstract 

, We give an explicit formula for the Hankel transform of a regular sequence in terms 

of the coefficients of the associated orthogonal polynomials and the sequence itself. 

' We apply this formula to some sequences of combinatorial interest, deriving interesting 

. combinatorial identities by this means. Further insight is also gained into the structure 

• ' of the Hankel transform of a sequence. 

^ ! 

d I Let /i„ denote a sequence of numbers, with /iq = 1. The Hankel transform [16] of this 

sequence is the sequence hn of determinants 

l/^i+J |o<i,jr'<n- 

00 ■ We shall assume that the sequence /i„ is regular (or "Catalan-like"), that is, /i„ 7^ for any 

^ . n. Much research [4, 6, 9, 10, 13, 19, 20, 21, 22] has focused on finding closed form expres- 

ly^ i sions for the Hankel transforms of important sequences. Some well-known sequences have 

^ I Hankel transforms with simple formulas (such as our examples below), but in general Hankel 

O ' transforms can have quite complex forms [9], or be totally resistant to current techniques of 

elucidation. The sequences we shall use in this note are susceptible to the following classical 
approach. If [14, 15] fin has a generating function g{x) expressible in the form 

X ■ f ^ /^O 

H ■ 9{x) - 



1 — aox 



\2X 

aix 



Xsx' 

a2X — 



— a^x — ■ ■ ■ 
then we have [14] 

n 

fc=i 

We can associate a sequence of monic orthogonal polynomials [5, 11, 25] Pn{x) with the 
sequence The /x„ then correspond to the moments of this family, for the linear functional 
C defined by 

C{x^) = flu. 
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The family {P„(x)} satisfies a three term recurrence relation 

Pn+iix) = (x - a„)P„(x) - A„P„_i(x), n>0, 

1. We can construct the P„(x) explicitly as follows [12]. 



where AoP_i(a:) = and Pi(x 
Let 



A,-, 



/ii 



/ii 

/i2 



DJx) 



/iO /il 
/il /i2 



1 X 



/ii+1 



and let 
Then 



P)„ = A. 



PJx) 



DJx) 



is the monic orthogonal sequence for C The orthogonality of the polynomials is equivalent 
to 

^{PmPn) = A1A2 ■ ■ ■ A„(5 



In particular, 

However, we know that [20] 
Thus we have that 



/:(p2) = AiA2---A„. 



hn — Y\ -^1-^2 ■ ■ ■ Aa 
fe=0 



hn = ll^{P',). 

k=0 

We wish now to evaluate C{P^). To do this, we let (a„^fc) be the matrix of coefficients of the 
polynomials {P„(x)}. Since the degree of P„(x) is n, this is a lower triangular matrix, with 
I's on the diagonal since the polynomials are monic. It is thus invertible (and if its elements 
are integers, so are those of the inverse), and we know that [1] the first column of its inverse 
is composed of the fin- Using the Cauchy formula for the coefficients of the product of two 
polynomials, along with the fact that £ operates on polynomials as follows: 



j=0 



j=0 



we obtain the following formula: 



2k 



^Pk) 



I x]«fcj«fc,i-i I /ii- 

i=0 \j=0 



(2) 
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Thus we obtain the following formula for the Hankel transform hn in terms of the elements 
of the coefficient array for the associated orthogonal polynomials: 

n 2k i 

^ n (^k,jak,i^j)fii. (3) 

k=0 1=0 j=0 

We state this proposition. 

Proposition 1. Let fin be a regular sequence, whose associated family of orthogonal polyno- 
mials {Pn{x)}n>o has coefficient array icLn,k)n k>o- Then the Hankel transform hn of fin is 
given by 

n 2k i 

^ TT y^('^ ak,jak,i-j)fJ'i- 

k=0 1=0 j=0 

Note also that the elements of this product, namely A1A2 ■ ■ ■ A^, are precisely the elements 
of of D where [17, 18, 26] the Hankel matrix H of fin has the LDU decomposition 

= {fi,+,) = LDL^, 

where L is the inverse of (ajj). Thus the elements of D are given by 

2n i 

^(^a„ja„,i_j)/ii, n = 0,1,2, 

We state this clS cl corollary to the proposition. 

Corollary 2. Let fin be a regular sequence, with Hankel matrix H = {fJ'i+j)i j^q- Then 

H = LDL^, 

where the elements dn of the diagonal matrix D are expressible as 

2k i 

dk = a'k,jO'k,i-j)fi'i- 

1 Examples 

In this section, we give some examples, concentrating on "simple" sequences of combinatorial 
significance whose Hankel transforms are well known. In each case the coefficient array of 
the associated family of monic orthogonal polynomials is a Riordan array [7, 23, 24]. These 
examples further explore links [1, 2, 3] between Riordan arrays, orthogonal polynomials and 
Hankel transforms. For those not familiar with Riordan arrays, we recall that the Riordan 
group [23, 24], is a set of infinite lower-triangular integer matrices, where each matrix is 
defined by a pair of generating functions g{x) = l+gix+g2x'^ + - ■ ■ and f{x) = fix+f2X^ + - ■ ■ 
where /i 7^ [24]. We assume in addition that /i = 1 in what follows. The associated matrix 



3 



is the matrix whose i-th column is generated by g{x)f{xy (the first column being indexed 
by 0). The matrix corresponding to the pair g, f is denoted by {g, f) or 7l{g, /). The group 
law is then given by 

{gJ)-{hJ) = {gJ){h,l) = {g{hof)Jof). 

The identity for this law is / = (l,x) and the inverse of {g,f) is {g,f)~^ = (1/(5' o Z),/) 
where / is the compositional inverse of / (defined by f{f{x)) = f{f{x)) = x)). 

Example 3. The binomial matrix (Pascal's triangle) which is the matrix with general term 
(2) is the Riordan array 

\ 1 — X 1 — X 

It has inverse B^-*^ = ((— given by 

1 X 



B 1 



X 1 — X 



The exponential Riordan group [7, 8], is a set of infinite lower-triangular integer matrices, 
where each matrix is defined by a pair of generating functions g{x) = go + gix + g2x'^ + ■ ■ ■ 
and f{x) = fiX + /2X^ + ■ ■ ■ where 7^ and /i 7^ 0. In what follows, we shall assume 

90 = fi = 1- 

The associated matrix is the matrix whose i-th column has exponential generating function 
g{x)f{xy/i\ (the first column being indexed by 0). The matrix corresponding to the pair 
/, g is denoted by [g, /]. The group law is given by 

[gJ]-[hJ] = [g{hof)Jof]. 

The identity for this law is / = [1, x] and the inverse of [g, f] is [g, f]^^ = [1/(5 o /), /] where 
/ is the compositional inverse of /. 

Example 4. The binomial matrix B = ( (^) ) is also an element of the exponential Riordan 
group. We have 

B = [e^x] and B^^ = [e-^x] . 

We are now is a position to apply our proposition to some sequences of combinatorial 
interest. 

Example 5. The Catalan numbers A000108 . It is well known that for the Catalan numbers 
Cn = ;^^(^^) we have hn = I for all n. This can be seen since the generating function 

g{x) = ^~^2x^^ Catalan numbers can be expressed as 

g{x) = i ^ . 

X 



The corresponding monic orthogonal polynomials 



fc=0 ^ ^ 



have coefficient array (^j^, (jq^j^j- We have 

P„+i(x) = (x-2)P„(x)-P„_i(x), Po{x) = l, Piix)=x-1. 
Equation (3) takes the form 

k=0 i=0 j=0 \ J / \ \ J I 



In fact in this case, since Aj = 1, we obtain the identity 

'^^ * ^k + j\fk + i— j 



i=0 j=Q 



Example 6. The central binomial numbers A000984 . The Hankel transform of the central 
binomial numbers (^^) is given by /i„ = 2". This can be seen since the generating function 
g(^x) = of the central binomial numbers can be expressed as 



2x' 

l-2x 



l-2x- 



2 

X 



1 - 2x 



x^ 



1 - 



The corresponding orthogonal polynomials Pj^x) have coefficient array i^ji^-, (jqi^^ • In this 
case, we have 

P„+i(a;) = (x - 2)P„(x) - P„_i(x), Po(x) = 1, P^{x) =x-2. 
The general term an,k of the Riordan array ^yT^' IJ+xf) given by 

( (n + k\ (n + k-l\\ , fn + k\ 2n + 0"+^ , 

Here we have used the notation 0" to denote the sequence with term 1, 0, 0, 0, . . . and gener- 
ating function 1. Thus 0"^ = 5o,n = [^^ = 0]. With this value for an,k, we then obtain 



n 2k i y„ .\ 



k=0 i=0 j=0 



Now Ai = 2, Aj = 1 for i > 1, and so we also have 

2k i 



(1 '\ 



j=0 i=0 

Thus we have, for instance, 

2k 

.^0 ^.^0 V 2 j y V - j) J k + j + 0''+^ k + j + 0''+'-^ V ^ 

Example 7. The large Schroder numbers A006318 . The large Schroder numbers Sn = 
^^=0 {^2k)^'' have Hankel transform /i„ = 2^ 2 ). This can be seen since the generating 
function g{x) of the large Schroder numbers can be expressed as 

1 

9{x) 



2x' 

l-2x 



2x 

l-3x 



2x2 

1 -3x 



1 

The corresponding monic orthogonal polynomials Pn{x) satisfy 

P„+i(x) = {x- 3)Pn{x) - 2P„„i(x), Po{x) = 1, Pi(x) = X - 2, 
with coefficient array equal to the Riordan array 

1 X 



For this array, we have 



l + 2x' l + 3x + 2x2 / ' 



(^n,k 

j=0 

With this value for a„,fc, we then have 

n 2k i 

/i„ = 2( 2 ) = Yl'^C^akjak,i-j)Si 

k=0 i=0 j=0 

and 

2k i 

^C^akjak,i-j)Si = 2^. 
Explicitly, this last equation gives us 

i=0 j=0 i=0 ^ ^ ^ m=0 \ -J ' ^ 
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Example 8. The factorial numbers A000142 . The factorial numbers n\ have Hankel trans- 
form hn = nr=o^'^- This can be seen since the generating function g{x) of the factorial 
numbers can be expressed as 



1 



9[X) 



1 — X — 



X 



1 — 3x 



1 — 5x 



9x' 



The corresponding monic orthogonal polynomials Pn{x) satisfy 

Pn+i(x) = (x - (2n + l))P„(x) - n^Pn^iix), Po{x) = 1, Pi(x) = X - 1, 
with coefficient array equal to the exponential Riordan array 

1 



X 



1 + X 1 + X 



with general term a„ ^ given by 



n\n\ 



k k\ 



\n—k 



With this value for a„ fc, we then obtain 



n 2k i n 
A;=0 i=0 j=0 i=0 



and 



2k i 
1=0 j=0 

Writing these expressions explicitly, we get 



and 



We deduce that 



n 2k i 
k=0 i=0 j=0 

2k i 
i=0 j=0 
2k i 

E(E(-i) 

i=0 j=0 



k\ f k 



2 J V- JJ J'-[^-J 



i=0 



k\ f k 



J J v-jJjK^-j)- 



k\f k 



)i\ = k\\ 



HI = 1, 



J/ V- J/ J-[^- J)- 
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or 



2k i 
i=Q j=0 



k 



1. 



Expanding the inner term in i and k as an array, we see that this expression represents the 
row sums of the array that begins 



-2 2 

-4 10 -12 6 

-6 24 -56 78 

-8 44 -152 346 

-10 70 -320 1010 



-60 20 
-520 500 
-2252 3560 



-280 70 
-3920 2870 



-1260 252 



Example 9. The derangement (or rencontres) numbers A000166 . It is well-known [16] that 
a sequence and its binomial transform, or its inverse binomial transform, have the same 
Hankel transform. Thus the derangement numbers [22] 



4=x:(-i)»-'(^)«, 

fc=0 ^ 



which are the inverse binomial transform of the factorial numbers n\, have Hankel transform 
hn = nr=o^'^- This may be seen from the continued fraction expansion of their generating 
function, given by 

1 



X 



1 - 2x 



4x^ 



Ax 



9x^ 



In this case, the coefficient array of the associated orthogonal polynomials is equal to that 
of the factorial numbers, multiplied on the left by the binomial matrix. We then obtain 



\e^x] 



X 



1 + x' 1 + 



X 



X 



X 1 



X 



as the Riordan array expression for the coefficient array of the associated orthogonal poly- 
nomials, with general term an,k given by 



Q"n,k 



E 



k kv ' 



With this value for an,fc, we then deduce, for instance, that 

2k i i 



i=0 j=0 



(=0 
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